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We  have  been  working  on  identifying  topoiogicai  phases  in  iridate  compounds  in  which  an  intricate  interpiay  between  spin-orbit 
interactions,  crystai  fieid  spiitting  and  couiomb  repuision  of  5d  eiectrons.  We  had  two  significant  achievements  under  this  grant. 
First,  we  used  time-resoived  opticai  spectroscopy,  which  is  highiy  sensitive  to  the  existence  of  energy  gaps,  to  study  the 
temperature  evoiution  of  the  eiectronic  structure  of  Sr2ir04  .  Taking  advantage  of  quaiitativeiy  distinct  reiaxation  dynamics  of 
photoexcited  carriers  exhibited  by  gapped  and  gapiess  systems,  we  find  a  ciear  change  in  the  uitrafast  dynamics  across  TN 
indicating  a  gap  opening  concomitant  with  antiferromagnetic  order.  Anaiysis  of  the  iong-time  reiaxation  dynamics  further  reveais 
behavior  that  is  consistent  with  in-gap  spectrai  weight  being  continuousiy  transferred  to  high  energies  through  TN  ,  supporting  a 
phase  transition  from  iow-temperature  insuiator  to  high-temperature  paramagnetic  metai  with  both  Mott-  Hubbard  and  Siater- 
type  characteristics. 

Secondiy,  we  observed  the  confinement  of  singie  particie  excitations  and  a  confinement  to  deconfinement  transition  in  a 
iayered,  quasi-2D,  strongiy  frustrated  antiferromagnetic  Na2lr03  iridate  system.  Recentiy  severai  research  groups  conjectured 
that  these  materiais  might  be  an  experimentai  reaiization  of  the  Kitaev  modei  of  a  spin  iiquid.  Aithough  it  was  experimentaiiy 
shown  that  this  system  does  undergo  ordering  at  reiativeiy  iow  temperatures  (T=15K)  and  hence  is  not  a  spin  iiquid,  it  was 
demonstrated  it  is  stiii  a  highiy  frustrated  spin  system. 

For  the  detaiis  of  these  achivements,  piease  refer  to  the  attached  papers.  We  are  aiso  inciuding  a  short  summary  of  the 
unpubiished  resuits  describing  work  done  after  our  iatest  pubiication. 
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Observation  of  a  metal-to-insulator  transition  with  both  Mott-Hubbard  and  Slater  characteristics 

in  Sr2lr04  from  time-resolved  photocarrier  dynamics 

D.  Hsieh,'  F.  Mahmood,*  D.  H.  Torchinsky,'  G.  Cao,^’^  and  N.  Gedik' 

^Department  of  Physics,  Massachusetts  Institute  of  Technology ,  Cambridge,  Massachusetts  02139,  USA 
^Center  for  Advanced  Materials,  University  of  Kentucky,  Lexington,  Kentucky  40506,  USA 
^Department  of  Physics  and  Astronomy,  University  of  Kentucky,  Lexington,  Kentucky  40506,  USA 
(Received  20  March  2012;  published  18  luly  2012) 

We  perform  a  time-resolved  optical  study  of  Sr2lr04  to  understand  the  influence  of  magnetic  ordering  on  the 
low  energy  electronic  structure  of  a  strongly  spin-orbit  coupled  J^-ff  —lf2  Mott  insulator.  By  studying  the  recovery 
dynamics  of  photoexcited  carriers,  we  find  that  upon  cooling  through  the  Neel  temperature  Tn  the  system  evolves 
continuously  from  a  metal-like  phase  with  fast  (~50  fs)  and  excitation  density  independent  relaxation  dynamics 
to  a  gapped  phase  characterized  by  slower  (~500  fs)  excitation  density-dependent  bimolecular  recombination 
dynamics,  which  is  a  hallmark  of  a  Slater-type  metal-to-insulator  transition.  However  our  data  indicate  that 
the  high  energy  reflectivity  associated  with  optical  transitions  into  the  unoccupied  /etf  =1/2  band  undergoes 
the  sharpest  upturn  at  Tn,  which  is  consistent  with  a  Mott-Hubbard  type  metal-to-insulator  transition  involving 
spectral  weight  transfer  into  an  upper  Hubbard  band.  These  findings  show  Sr2lr04  to  be  a  unique  system  in 
which  Slater-  and  Mott-Hubbard-type  behaviors  coexist  and  naturally  explain  the  absence  of  anomalies  at  Tn  in 

BAGS  number(s):  71.27. -|-a,  71.70.Ej,  78.47. -p 

exactly  at  7/v.'^  Therefore  the  relevant  experimental  question  is 
whether  heating  above  7/v  brings  Sr2lr04  into  a  paramagnetic 
insulating  phase  or  into  a  paramagnetic  metallic  phase  and 
what  the  order  of  the  MIT  is. 

Owing  to  an  absence  of  clear  anomalies  at  in 
transport, thermodynamic,'^  and  optical  conductivity 
data,'®  there  have  been  conflicting  interpretations  about  how 
the  insulating  gap  behaves  across  Tn-  In  this  article  we  use 
time-resolved  optical  spectroscopy,  which  is  highly  sensitive 
to  the  existence  of  energy  gaps,^"’^'  to  study  the  temperature 
evolution  of  the  electronic  structure  of  Sr2lr04.  Taking 
advantage  of  qualitatively  distinct  relaxation  dynamics  of 
photoexcited  carriers  exhibited  by  gapped  and  gapless  systems, 
we  And  a  clear  change  in  the  ultrafast  dynamics  across  Tn 
indicating  a  gap  opening  concomitant  with  antiferromagnetic 
order.  Analysis  of  the  long-time  relaxation  dynamics  further 
reveals  behavior  that  is  consistent  with  in-gap  spectral  weight 
being  continuously  transferred  to  high  energies  through  Tn, 
supporting  a  phase  transition  from  low-temperature  insulator 
to  high-temperature  paramagnetic  metal  with  both  Mott- 
Hubbard  and  Slater-type  characteristics. 

II.  EXPERIMENTAL  DETAILS 

In  our  experiment  we  used  a  Ti:  sapphire  oscillator  pro¬ 
ducing  laser  pulses  with  center  wavelength  795  nm  {hv  — 
1.56  eV)  and  near  80  fs  duration.  The  fluence  of  the  pump 
pulse,  which  excites  electrons  across  the  insulating  gap 
from  the  occupied  7eff  =  3/2  and  Teff  =1/2  bands  to  the 
unoccupied  /eft  =1/2  band  [inset  Fig.  1(a)],  was  varied  with 
neutral  density  filters  to  tune  the  photoexcited  carrier  density 
while  the  probe  fluence  was  maintained  at  4  /xJ/cm^.  Both 
beams  were  focused  to  a  70-/xm  full  width  at  half  maximum 
(FWHM)  spot  on  the  (001)  cleaved  face  of  the  sample.  The 
80-MHz  repetition  rate  was  reduced  to  1.6  MHz  with  a  pulse 
picker  to  eliminate  steady-state  heating  of  the  sample.  Use  of  a 


transport  and  thermodynamic  measurements. 

DOT  10.1103/PhysRevB. 86.035128 

I.  INTRODUCTION 

Iridium  oxides  are  unique  5cl  electronic  systems  in  which 
spin-orbit  coupling,  electronic  bandwidth  (W),  and  on-site 
Coulomb  interactions  (U)  occur  on  comparable  energy  scales. 
Their  interplay  can  stabilize  a  novel  /eft  =1/2  Mott  insulating 
state  in  which  a  correlation  gap  is  opened  by  only  moderate 
Coulomb  interactions  owing  to  a  spin-orbit  coupling  induced 
band  narrowing.*  Depending  on  the  underlying  lattice,  this  in¬ 
sulating  state  is  predicted  to  realize  a  variety  of  exotic  quantum 
phases  including  antiferromagnetic  =1/2  Mott  insulators 
on  the  perovskite  lattice,*"^  correlated  topological  insulators 
and  semimetals  on  the  pyrochlore  lattice,"*"®  and  topological 
spin  liquids  on  the  hyper-kagome^  and  honeycomb  lattices.* 
Even  more  tantalizing  possibilities  are  predicted  to  occur 
upon  chemically  doping  these  systems,  ranging  from  high-!/ 
superconductivity®’"*  to  spin-triplet  superconductivity.* 

Intensive  research  has  been  conducted  on  the  layered 
perovskite  iridate  Sr2lr04  owing  to  its  structural  and  electronic 
similarities  to  undoped  high- 7/  cuprates  such  as  La2Cu04.*’® 
The  ground-state  electronic  structure  of  Sr2lr04  consists 
of  a  completely  filled  band  with  total  angular  momentum 
=  3/2  and  a  narrow  half-filled  Teff  =1/2  band  near  the 
Fermi  level  The  latter  is  split  into  an  upper  Hubbard 

band  (UHB)  and  lower  Hubbard  band  (LHB)  due  to  on-site 
Coulomb  interactions  and  exhibits  antiferromagnetic  ordering 
of  the  effective  Jeff  =1/2  moments  below  a  Neel  temper¬ 
ature  Tn  =  240  K*  analogous  to  La2Cu04.  Although  this 
insulating  ground  state  has  been  established  by  angle-resolved 
photoemission  spectroscopy'  and  resonant  x-ray  scattering" 
measurements,  whether  Sr2lr04  is  a  Mott-type  (U  )$>  W) 
insulating  phase  typical  of  3d  transition  metal  oxides  or  a 
Slater-type  (U  ^  W)  insulating  phase  is  experimentally  un¬ 
known  and  remains  a  topic  of  active  theoretical  debate."*’**’*"* 
Whereas  a  Mott-Hubbard-type  metal-to-insulator  transition 
(MIT)  is  discontinuous  and  occurs  at  temperatures  greater 
or  equal  to  Tn,  a  Slater-type  MIT  is  continuous  and  occurs 
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FIG.  1.  (Color)  Normalized  time-resolved  reflectivity  traces 
(AR/R)norm  of  Sr2lr04  collected  in  the  temperature  ranges  (a)  4- 
120  K,  (b)  120-210  K,  (c)  210-300  K,  and  (d)  300-350  K.  Curves  are 
collected  in  1 0-K  intervals  with  a  pump  fluence  of  1 5 .4  /r  J/cm^ .  Inset 
of  panel  (a)  shows  a  schematic  of  the  low-energy  electronic  density 
of  states  (DOS)  based  on  calculations.'^  The  red  arrow  denotes  the 
optical  transition  being  excited  by  the  pump  pulse.  Inset  of  panel 
(d)  shows  the  T  =  300  K  trace  out  to  100-ps  time  delay. 

double-modulation  scheme^®  provided  sensitivity  to  the  frac¬ 
tional  change  of  reflectivity  on  the  order  AR/R  ~  10“^.  Single 
crystals  of  Sr2lr04  were  grown  using  a  self-flux  technique 
and  magnetization  measurements  show  a  magnetic  ordering 
temperature  at  7jv  =  240 

III.  RESULTS  AND  ANALYSIS 

Figure  1  shows  typical  time-resolved  reflectivity  transients 
measured  over  a  range  of  temperatures  spanning  4-350  K, 
which  have  all  been  normalized  to  their  negative  peak  values 
in  order  to  emphasize  the  recovery  dynamics.  We  note  that 
the  temperature  here  refers  to  that  before  the  pump  excitation, 
which  will  be  lower  than  the  instantaneous  electronic  temper¬ 
atures  reached  in  the  ps  time  windows  following  photoexcita¬ 
tion.  Following  the  pump  excitation,  the  electronic  temperature 
is  elevated  above  the  initial  equilibrium  temperature  and  all 
AR/R  traces  exhibit  a  rapid  negative  spike,  which  indicates 
a  decrease  in  reflectivity.  Within  approximately  1  ps,  the 
reflectivity  recovers  to  a  small  negative  offset  that  persists 
beyond  100  ps  [Fig.  1(d)  inset].  By  performing  a  temperature 
dependence  study  we  And  that  the  initial  ~l-ps  recovery 
can  actually  be  decomposed  into  two  separate  components 
as  follows.  Figure  1(a)  shows  that  between  4  and  120  K 
there  is  no  discernible  change  in  the  recovery  dynamics.  Upon 
heating  from  120  to  210  K  [Fig.  1(b)],  we  find  no  temperature 
dependence  within  the  first  ~  100  fs  but  observe  a  clear  slowing 
down  of  recovery  dynamics  after  this  time.  Further  heating  of 
the  sample  from  210  to  300  K  again  causes  no  change  in  the 
first  recovery  component  but  causes  the  second  component  to 
become  faster  [Fig.  1(c)].  All  temperature  dependence  then 
completely  shuts  off  above  300  K  [Fig.  1(d)].  To  understand 
the  physical  processes  underlying  these  trends,  we  perform 
detailed  fits  to  the  AR/R  traces. 

We  fit  the  un-normalized  time-resolved  reflectivity  tran¬ 
sients  [A/?//?](t)  at  all  temperatures  to  a  convolution 


[/  <S>  g](0  =  fZo  of  a  bi-exponential  decay 

function  f{t)  and  a  Gaussian  instrument  resolution  function 
git)  where 


/(0  = 


0;  if  f  <  0 

Ae-'/^'  -b  +  C;  if  f  >  0, 


and 


git)  = 


1 

ay/ In 


where  a  is  experimentally  determined  by  measuring  an 
autocorrelation  of  our  laser  pulses. 

To  demonstrate  that  the  [Ai?//?](t)  traces  are  described  by 
a  minimum  of  two  exponential  decays,  we  fit  [/  ®  g](f)  to 
our  data  at  positive  time  delays  to  both  a  single  and  double 
exponential  decay.  We  show  fits  to  the  data  in  two  extreme 
cases:  (i)  data  taken  at  a  temperature  T  —  300  K  and  pump 
fluence  F  —  15.4  /zJ/cm^,  where  the  magnitude  of  the  slower 
decay  component  is  shown  to  be  smallest  [shown  later  in 
Fig.  4(c)],  (ii)  data  taken  at  a  temperature  T  —  180  K  and 
pump  fluence  F  —  9  /xJ/cm^,  where  the  magnitude  of  the 

slower  decay  component  is  shown  to  be  largest  [shown  later  in 
Fig.  4(c)].  Figure  2  shows  that  in  both  extremes,  the  data  can 
only  be  satisfactorily  described  by  two  exponential  decays. 
A  decomposition  of  the  fits  into  its  individual  components 
(Fig.  2)  shows  that  the  slow  component  is  required  to  account 
for  the  curvature  in  AR/R  around  0.5  ps.  The  fitted  parameters 
listed  in  red  show  that  their  standard  error  is  small,  less  than 
the  size  of  the  symbols  plotted  in  all  the  figures.  We  found  that 
the  fit  results  are  not  sensitive  to  changes  of  order  ±40%  in 
their  initial  estimated  values. 
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FIG.  2.  (Color)  Transient  reflectivity  traces  taken  at  T  =  300 
K  and  F  =  15.4  /xj/cm^  and  fitted  to  a  (a)  single  and  (b)  double 
exponential  decay.  Transient  reflectivity  traces  taken  at  T  =  180 
K  and  F  =  1.9  /xj/cm^  and  fitted  to  a  (c)  single  and  (d)  double 
exponential  decay.  Curves  show  the  overall  fit  as  well  as  the  separate 
components  of  the  fit.  The  fitted  parameters  and  their  associated  fitting 
errors  (in  parentheses)  are  listed  in  red. 
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FIG.  3.  (Color)  (a)  A  typical  un-normalized  time-resolved  reflec¬ 
tivity  transient  (AR/R)  of  Sr2lr04  overlaid  with  a  bi-exponential  fit 
described  in  the  text  (red  curve).  Inset  shows  the  Gaussian  instrument 
time  resolution  (red  curve)  superimposed  on  300-K  data  where  decay 
dynamics  are  dominated  by  the  fast  (ti)  component  [Fig.  4(b)].  This 
shows  that  the  fast  initial  decay  is  not  resolution  limited,  (b)  The  fitted 
amplitude  (A)  of  the  fast  exponential  decay  component  measured  over 
a  range  of  pump  fluences  corresponding  to  to  ~10“^  pump 

photons  per  iridium  atom.  The  red  line  is  a  linear  fit.  (c)  Temperature 
dependence  of  the  fitted  decay  time  (ti)  of  the  fast  exponential  decay 
component  measured  over  a  range  of  pump  fluences.  (d)  Schematic 
showing  that  the  initial  fast  decay  process  is  governed  by  energy 
relaxation  of  photoexcited  electrons  and  holes  towards  the  band  edges 
via  optical  phonon  emission.  The  errors  in  the  fit  parameters  are 
smaller  than  the  size  of  the  symbols. 


A.  Fast  decay  component 

We  begin  by  trying  to  understand  the  physical  origin  of 
the  first  decay  component,  which  is  resolved  by  our  ultrashort 
laser  pulses  [Fig.  3(a)].  To  investigate  whether  this  fast  decay 
arises  from  thermalization  of  photoexcited  carriers  via  carrier- 
carrier  scattering,  whose  rate  should  depend  on  the  number  of 
photoexcited  carriers, we  measured  ti  over  a  range  of  pump 
fluences  spanning  a  photoexcitation  density  d  between  ~10“^ 
to  ~10“^  photons  per  iridium  site.  We  estimate  d  using  the 
expression  d  =  where  N  is  the  number  of  photons  per 
pump  pulse,  n  is  the  number  of  iridium  sites  per  unit  volume, 
a  is  the  area  of  the  sample  illuminated  by  the  pump  pulse,  and 
f  is  the  optical  penetration  depth  of  photons  with  wavelength 
k  — 195  nm.  The  penetration  depth  §  is  derived  through  the 
Beer-Lambert  law  f  The  extinction  coefficient  k, 

equal  to  the  imaginary  part  of  the  index  of  refraction,  can 
be  calculated  using  Eq.  (8)  in  Ref.  24  from  the  amplitude  (R) 
and  phase  (0)  of  the  reflection  coefficient  of  Sr2lr04.  We  use 
available  published  values  of  R  and  </>  from  lightly  Rh-doped 
Sr2lr04  thin  films, whose  optical  conductivity  spectrum  is 
nearly  identical  to  bulk  single  crystalline  Sr2lr04. 

A  linear  behavior  of  the  component  amplitude  (A)  over 
this  fluence  range  [Fig.  3(b)]  shows  that  the  number  of  pho¬ 
toexcited  carriers  is  indeed  proportional  to  the  pump  fluence. 
Figure  3(c)  shows  that  ti  ss  50  fs  exhibits  no  discernible 
fluence  nor  femperature  dependence  between  4  and  350  K. 


This  implies  that  the  initial  decay  is  not  caused  by  photocarrier 
thermalization  but  more  likely  by  photocarrier  cooling,  and 
that  this  cooling  is  not  mediated  by  thermally  occupied 
phonons.  The  participation  of  magnons  is  also  negligible 
because  ti  is  insensitive  to  7^.  Given  that  the  fast  time 
scale  of  ti  is  consistent  with  typical  optical  phonon-mediated 
cooling  processes  and  that  the  Debye  temperature  of  Sr2lr04'® 
far  exceeds  our  measurement  temperatures,  we  conclude  that 
the  initial  fast  recovery  of  AR/R  is  due  to  the  cooling  of 
photoexcited  carriers  via  generation  of  hot  optical  phonons 
[Fig.  3(d)]. 


B.  Slow  decay  component 

The  initial  fast  recovery  of  AR/R  is  followed  by  a  second 
slower  decay  component,  whose  amplitude  (B)  exhibits  an 
upturn  upon  cooling  the  starting  equilibrium  temperature  of 
the  system  through  7jv  and  then  ceases  to  grow  further  below 
a  temperature  Tq  =  175  K  [Fig.  4(a)].  The  amplitude  of  the 
fast  component  (A)  shows  comparatively  little  temperature 
dependence  [Fig.  4(a)],  therefore  the  temperature  dependence 
of  their  relative  amplitude  (B /A)  is  dominated  by  the  features 
observed  in  B  [Figs.  4(b)  and  4(c)].  This  indicates  that  at 
temperatures  far  above  7)v  the  electronic  energy  relaxation 
occurs  predominantly  through  the  lone  process  of  optical 
phonon  generation  described  in  Fig.  3  and  that  a  separate 
relaxation  mechanism  grows  near  7)v  before  saturating  below 
Tq.  The  decay  time  r2  shows  a  similar  upturn  upon  cooling 
through  7]v  from  around  0.36  ps  at  350  K  to  0.55  ps  at  Tq 
and  then  exhibits  a  marked  change  in  temperature  dependence 
below  To  [Fig.  4(d)].  A  sharp  rise  in  relaxation  time  at  7)v 
typically  signifies  the  development  of  an  energy  gap,^°7i 
where  a  depletion  of  states  around  Ef  greatly  reduces  the 
efficiency  of  photocarrier  relaxation  from  above  to  below  E  f . 

The  hallmark  of  a  fully  formed  energy  gap  is  that  photoex¬ 
cited  occupied  states  above  Ef  and  empty  states  below  Ef 
can  only  combine  in  a  pairwise  fashion.  Therefore,  unlike 
in  a  metal,  the  recombination  rate  should  be  proportional 
to  the  density  of  photoexcited  electron-hole  pairs. We  find 
that  a  clear  fluence  dependence  of  r2  develops  only  below 
To  [Fig.  4(d)]  and  that  the  low-temperature  relaxation  rate 
Tj"'  indeed  increases  linearly  with  fluence  below  5  /rJ/cm^ 
[Fig.  4(e)].  This  shows  that  below  Tc,  t2  represents  the 
time  scale  for  recombination  of  photoexcited  electron-hole 
pairs  across  the  insulating  gap,  which  in  the  Mott-Hubbard 
limit  corresponds  to  the  recombination  of  empty  and  dou¬ 
bly  occupied  sites  [inset  Fig.  4(e)].  This  ultrafast  sub-ps 
recombination  time  in  Sr2lr04  is  comparable  to  that  observed 
in  the  two-dimensional  (2D)  antiferromagnetic  cuprate  Mott 
insulators,^^  which  has  been  attributed  to  additional  recombi¬ 
nation  channels  involving  hot  magnon  generation.  Altogether 
these  observations  allow  us  to  conclude  that  at  temperatures  far 
above  7)v  the  system  behaves  like  a  metal.  Near  Tn  the  density 
of  states  around  Ff  is  continually  depleted  upon  cooling, 
which  is  consistent  with  optical  conductivity  measurements,'^ 
and  the  insulating  gap  is  fully  developed  below  Tq  . 

To  rule  out  the  possibility  that  the  observed  temperature- 
dependent  features  are  artifacts  of  cross  correlation  between 
our  five  fitting  parameters,  we  note  that  based  on  Figs.  3(c)  and 
4(a)  the  parameters  A  and  Ti  exhibit  no  pronounced  features 
at  the  temperatures  of  interest  7)v  and  Tq.  Moreover,  we  show 
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FIG.  4.  (Color)  (a)  Temperature  dependence  of  the  fit  parameters 
A  and  B,  normalized  to  their  values  at  T  —  350  K,  taken  in  the 
low  pump  fluence  regime,  (h)  Temperature  dependence  of  B/A 
derived  from  data  shown  in  (a),  (c)  Temperature  dependence  of  the 
amplitude  ratio  and  (d)  the  decay  time  of  the  slower  component  (T2) 
measured  using  three  different  fluences  15.4  /rj/cm^,  1.9  /rj/cm^, 
and  0.18  /r  j/cm^.  (e)  Fluence  dependence  of  the  slower  decay  rate 
measured  at  4  K.  Straight  line  is  a  guide  to  the  eye  showing  a  linear 
dependence  at  low  fluences.  The  errors  in  the  fit  parameters  are 
less  than  the  size  of  the  symbols.  A  schematic  of  the  himolecular 
relaxation  process  involving  the  annihilation  of  photoexcited  empty 
and  doubly  occupied  sites  via  emission  of  optical  phonons  or  magnons 
is  shown  in  the  inset,  which  emphasizes  the  Mott-Hubbard  character 
of  the  system,  (f)  Schematic  showing  the  temperature  evolution  of 
the  low-energy  electronic  density  of  states  (DOS)  in  Sr2lr04  based 
on  our  data.  Filled  and  empty  circles  denote  the  photoexcited  electron 
and  hole  respectively. 


in  the  following  section  that  the  feature  at  Tn  observed  in  the 
temperature  dependence  of  C  is  completely  independent  of  the 
fitted  parameters  associated  with  the  two  exponential  decays. 
Therefore  the  only  parameters  that  show  features  atTui  and  Tg 
and  may  be  interdependent  are  B  and  12.  However,  the  fitting 
errors  associated  with  B  and  T2  are  very  small  (Fig.  2),  which 
is  further  support  that  their  temperature  dependence  is  real  and 
robust. 

The  occurrence  of  the  metal-to-insulator  transition  at  Tjv 
in  Sr2lr04  (Fig.  4)  distinguishes  it  from  archetypal  Mott 
insulators  such  as  MnO  where  the  insulating  gap  persists 
even  in  the  absence  of  long-range  magnetic  order.^^  The 


continuous  nature  of  the  MIT  and  its  development  over  a  broad 
(0.7  ^T/Tff  <  1.4)  temperature  window  further  precludes 
a  Mott-Hubbard  description,  which  predicts  sharp  first-order 
MITs  like  in  V203.^*  Although  our  measurements  do  not  rule 
out  gaps  beginning  to  form  in  microscopically  segregated 
regions  above  a  bulk  MIT  temperature  (7mit),  such  as  the 
case  in  doped  ¥203,^®  the  fact  that  such  gap  opening  takes 
place  exactly  across  in  Sr2lr04  defies  a  strict  Mott-Hubbard 
description.  We  rule  out  the  possibility  of  a  disorder  broadened 
Tmit  in  our  samples  based  on  their  sharp  magnetic  susceptibil¬ 
ity  curves. Although  these  results  point  towards  a  second- 
order  Slater-type  MIT,  Sr2lr04  does  not  conform  to  a  weakly 
correlated  (!/<$(  W)  spin-density  wave  description,  which  are 
driven  purely  by  Fermi  surface  nesting  and  are  typically  only 
partially  gapped  below  7^.  Unlike  a  conventional  spin-density 
wave  system,  magnetic  ordering  in  Sr2lr04  does  not  change 
the  size  of  the  unit  cell^'  and  local  moment  fluctuations  exist 
well  above  Tn  according  to  both  magnetic  susceptibility^’'^ 
and  magnetic  diffuse  x-ray  scattering  measurements.^"  While 
these  magnetic  signatures  provide  evidence  for  correlation 
physics  at  play,  a  direct  electronic  distinction  between  a 
correlation  driven  MIT  and  a  Fermi  surface  nesting  driven  MIT 
is  whether  or  not,  respectively,  spectral  weight  is  transferred 
from  low  (near  Ef)  energies  to  high  (of  order  U)  energies 
upon  traversing  the  MIT."' 

C.  Offset  component 

To  determine  whether  the  in-gap  spectral  weight  lost  at 
Tn  is  being  transferred  to  upper  and  lower  Hubbard  bands, 
we  investigate  the  temperature  dependence  of  the  optical 
reflectivity  at  the  energy  scale  of  optical  transitions  from  the 
occupied  /eff  =1/2  and  /eff  =  3/2  bands  to  the  unoccupied 
7gjj=l/2  band  [inset  Fig.  1(a)]  by  studying  the  small 
negative  offset  term  C  in  the  un-normalized  AR/R  traces.  The 
two  exponential  relaxation  processes  have  been  identified  as 
intraband  cooling  and  recombination  processes  that  both  lead 
to  a  transfer  of  energy  from  the  electronic  to  lattice  subsystems, 
which  brings  them  into  thermal  equilibrium  within  the  first  few 
ps  [Fig.  5(a)].  Subsequent  cooling  of  this  heated  sample  spot 
back  to  the  Initial  temperature  takes  place  through  the  diffusion 


FIG.  5.  (Color)  (a)  Diagram  illustrating  the  temporal  equilibra¬ 
tion  between  the  electronic  and  lattice  subsystems  following  a  pump 
pulse.  The  slowest  process  (C)  is  the  cooling  of  the  equilibrated 
system  via  diffusion  of  hot  carriers  or  hot  phonons  away  from  the 
measured  sample  spot,  (b)  Temperature  dependence  of  the  offset  (C) 
measured  using  a  pump  fluence  of  15.4  /rj/cm".  A  clear  minimum 
in  C  is  observed  at  Tn  corresponding  to  a  rapid  rise  in  reflectivity  at 
1.5  eV  at  Tn  (see  inset  schematic).  The  errors  in  the  fit  parameters 
are  less  than  the  size  of  the  symbols. 
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FIG.  6.  (Color)  Temperature  dependence  of  un-normalized  AR/R  traces  that  are  plotted  over  a  range  to  emphasize  the  long  time  offset. 
Data  are  plotted  from  (a)  T  —  4-130  K,  (b)  140-240  K,  and  (c)  250-350  K  in  increments  of  10  K.  (d)  The  value  of  AR/R  at  a  time  delay  of 
5  ps  [red  aiTows  in  panels  (a)-(c)]  plotted  as  a  function  of  temperature. 


of  hot  carriers  or  hot  phonons  away  from  the  laser  illuminated 
area,  which  is  estimated  to  well  exceed  100  ps  based  on  thermal 
conductivity  and  heat  capacity  data  for  Sr2lr04.^®  To  estimate 
the  time  required  for  heat  to  escape  the  laser  excited  region 
of  the  Sr2lr04  crystal,  we  calculate  its  thermal  diffusivity 
Dth  =  ^  where  k  is  the  thermal  conductivity,  p  is  the  density, 
and  c  is  the  specihc  heat.  We  use  the  values  measured  at 
r  =  100  K  in  Ref.  16  of  A:  =  2.5  W/m  K,  p  =  7440  kg/m^  and 
c  =  30  J/kg  K,  which  yields  Z)th  =  1.1  x  10“^  m^/s.  From  this 
we  estimate  the  time  needed  for  heat  to  diffuse  out  of  a  l-/rm^ 
area  to  be  around  1  ns.  Therefore  to  a  very  good  approximation 
the  offset  C  is  a  measure  of  the  fractional  change  in  reflectivity 
at  1.5  eV  due  to  a  small  temperature  change  AT  [Fig.  5(a)], 
namely  C  =  We  estimate  that  at  high  pump 

fluences  AT  is  of  order  1  K  because  C  is  of  order  10“"^  and 
the  fractional  decrease  in  optical  conductivity  at  1.5  eV  is 
estimated  to  be  2%  from  10  to  500  The  ability  to  resolve 
such  small  temperature-induced  changes  in  reflectivity  is  an 
advantage  our  nonequilibrium  technique  has  over  conventional 
equilibrium  optical  spectroscopy.  Figure  5(b)  shows  that  the 
temperature  dependence  of  C  exhibits  a  broad  negative  peak 
with  an  extremum  exactly  at  Tn,  which  indicates  that  a  small 
temperature  rise  causes  the  largest  decrease  of  the  1.5-eV 
reflectivity  exactly  at  7)v.  From  this  we  can  infer  a  broad  tem¬ 
perature  window  (0.6  <  T/T^  <  1.4)  within  which/?  (1.5  eV) 
starts  to  rise  sharply  with  cooling,  with  the  greatest  slope 
occurring  at  7jv  [inset  Fig.  5(b)].  The  fact  that  this  temperature 
window  largely  coincides  with  that  over  which  T2  increases 
most  drastically  [Fig.  4(d)]  is  consistent  with  a  Mott-Hubbard 
mechanism  where  spectral  weight  is  transferred  from  low  in¬ 
gap  energies  (~0.1  eV^^)  to  energies  (1.5  eV)  far  exceeding  it. 

To  demonstrate  that  this  behavior  is  independent  of  our 
fitting  procedure,  we  show  that  the  temperature  dependence 


of  the  long  time  offset  in  AR/R,  the  asymptotic  value  of 
the  un-normalized  reflectivity  variation  after  the  first  two 
fast  decay  processes,  can  be  obtained  even  without  fitting. 
Figure  6(a)  shows  that  between  T  —  A  and  130  K,  this 
asymptotic  value  does  not  show  a  measurable  systematic 
variation  with  temperature.  Upon  heating  from  140  to  240  K 
[Fig.  6(b)],  the  offset  clearly  decreases  with  heating  and  upon 
heating  further  from  250  to  350  K  [Fig.  6(c)],  the  offset 
exhibits  a  clear  increase  with  heating.  By  plotting  the  value 
of  AR/R  at  5  ps  as  a  function  of  temperature  [Fig.  6(d)], 
we  can  already  visualize  the  temperature  dependence  of  the 
offset,  which  is  very  similar  to  that  obtained  through  fitting 
[Fig.  5(b)].  This  proves  that  the  peak  observed  at  in  the 
offset  C  is  independent  of  how  we  fit  the  initial  fast  dynamics. 

We  note  that  variations  in  reflectivity  can  also  be  caused 
by  a  narrowing  or  shift  of  Lorentz  oscillators  representing 
direct  optical  transitions  at  1 .5  eV  that  do  not  involve  spectral 
weight  redistribution,  and  that  our  measurements  do  not 
distinguish  between  density  of  states  changes  in  the  occupied 
=  3/2  and  Jeff  =1/2  initial  states  versus  the  unoccupied 
=1/2  final  state.  However,  there  is  no  current  explanation 
for  how  such  changes  can  be  induced  via  magnetic  ordering 
in  the  7eff  =1/2  band.  Regardless,  the  observed  temperature 
dependence  of  the  quasiequilibrium  reflectivity  at  1.5  eV 
is  unusual  because  the  relaxation  rate  behavior  suggests  a 
Slater-type  behavior  (Fig.  4),  which  should  be  dominated  by 
low-energy  physics  near  the  Fermi  level  and  have  little  effect 
on  the  high-energy  reflectivity. 

IV.  CONCLUSIONS 

Our  results  taken  altogether  suggest  that  unlike  La2Cu04, 
which  remains  a  Mott  insulator  far  above  its  magnetic  ordering 
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temperature,  Sr2lr04  undergoes  a  metal-to-insulator  transition 
across  Tjy.  This  suggests  that  Sr2lr04  is  more  accurately 
described  as  an  intermediate  coupling  ([/  ss  W)  insulator  that 
exhibits  both  Mott-Hubbard  and  Slater-type  characteristics. 
Our  observation  that  this  metal-to-insulator  transition  takes 
place  over  a  wide  temperature  window  compared  to  7jv 
naturally  explains  the  lack  of  sharp  anomalies  at  Tn  in 
transport,  thermodynamic,  and  optical  conductivity  data.  We 
argue  that  these  may  be  signatures  of  a  rare  example  of  a 
temperature-controlled  continuous  metal-to-insulator  transi¬ 
tion  in  a  quasi-two-dimensional  system  hitherto  unobserved 
in  any  J-electron  material,^®  and  may  also  be  applicable  to 
the  wider  class  of  J^s  =1/2  Mott  insulating  perovskite, 
honeycomb,  hyperkagome,  and  pyrochlore  iridates.  Moreover, 
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transport  data  on  5d  Os  oxides  suggestive  of  continuous 
metal-to-insulators  transitions^^  show  that  the  confluence  of 
strong  spin-orbit  coupling  and  on-site  Coulomb  interactions 
is  a  general  playground  for  unconventional  metal-to-insulator 
transitions. 
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We  use  ultrafast  optical  spectroscopy  to  observe  binding  of  charged  single-particle  excitations  (SE) 
in  the  magnetically  frustrated  Mott  insulator  Na2lr03.  Above  the  antiferromagnetic  ordering  temperature 
(T ^)  the  system  response  is  due  to  both  Hubbard  excitons  (HE)  and  their  constituent  unpaired  SE.  The 
SE  response  becomes  strongly  suppressed  immediately  below  T We  argue  that  this  increase  in  binding 
energy  is  due  to  a  unique  interplay  between  the  frustrated  Kitaev  and  the  weak  Heisenberg-type  ordering 
term  in  the  Hamiltonian,  mediating  an  effective  interaction  between  the  spin-singlet  SE.  This  interaction 
grows  with  distance  causing  the  SE  to  become  trapped  in  the  HE,  similar  to  quark  confinement  inside 
hadrons.  This  binding  of  charged  particles,  induced  by  magnetic  ordering,  is  a  result  of  a  confinement- 
deconfinement  transition  of  spin  excitations.  This  observation  provides  evidence  for  spin  liquid  type 
behavior  which  is  expected  in  Na2lr03. 
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Spin-orbit  coupling  (SOC)  can  give  rise  to  bigbly 
nontrivial  physics.  Prime  examples  of  tbe  role  of  SOC 
in  condensed  matter  systems  are  tbe  topological  insulators 
wbicb  have  a  nontrivial  topology  of  tbeir  band  structure 
due  to  sufficiently  strong  SOC  [1,2].  In  tbe  case  of  simpler 
“band  topological  insulators,”  tbe  gap  is  determined  by 
tbe  spin-orbit  coupling  and  tbe  system  can  be  treated  as 
noninteracting.  These  materials  have  been  subject  to 
intensive  research  and  are  relatively  well  understood. 
Much  less  clear  is  the  situation  in  which  the  insulating 
state  before  the  “turning  on”  of  SOC  was  not  a  trivial 
band  insulator  but  an  insulator  with  a  gap  driven  by 
electton-electron  interactions,  such  as  a  Mott  insulator. 

Iridate  compounds  belong  to  the  class  of  materials  in 
which  electron-electron  interactions  play  an  essential  role. 
Expected  to  be  metallic  based  on  a  simple  electron  count, 
these  systems  are  insulators  exhibiting  Mott-type  behavior. 
On  the  other  hand,  due  to  the  extended  nature  of  the  5d 
orbitals,  the  on-site  Coulomb  repulsion  has  a  moderate 
value  (U  «  0.4-2. 5  eV)  and  SOC  («0.4-l  eV)  effectively 
competes  with  electron-electron  interactions  [3]. 

One  of  the  most  intriguing  proposals  of  novel  physics 
made  for  iridate  compounds  was  that  the  interplay  between 
spin-orbit  interactions,  crystal  field  splitting,  and  Coulomb 
repulsion  of  5d  electrons  in  Na2lr03  can  lead  to  a  formation 
of  effective  moments  with  =1/2  on  every  Ir-0 
octahedron  with  highly  anisotropic  nearest  neighbor  cou¬ 
pling.  This  coupling  has  a  very  special  form  and,  given  its 
layered  quasi-2D  honeycomb  lattice  structure,  was  pro¬ 
posed  [4-7]  to  be  a  solid-state  realization  of  the  Kitaev 
model  [8]  of  a  spin  liquid. 

The  real  Hamiltonian  of  Na2lr03  is,  however,  not  a  pure 
Kitaev  model  and  should  also  have  conventional  terms  such 


as  Heisenberg-type  exchange  interaction  between  effective 
moments.  Such  terms  generally  spoil  the  symmetry  of  the 
pure  Kitaev  model  and  result  in  an  ordered  ground  state 
[5,7,9]  in  the  limit  of  low  temperatures.  The  structure  of  the 
ground  state  should  then  depend  on  the  details  of  the  extra 
term  in  the  Hamiltonian.  Neutron  studies  [10-12]  have 
revealed  that  Na2lr03  has  an  antiferromagnetic  ground 
state  of  “zigzag”  type  [Eig.  3(b)]  with  a  Neel  temperature  of 
Tj^  =  15.3  K.  The  minimal  Hamiltonian  within  the  frame¬ 
work  of  the  modified  Kitaev  model  that  can  give  such  a 
ground  state  consists  of  an  antiferromagnetic  Kitaev  term 
and  a  ferromagnetic  Heisenberg  term  [5,7,9,13].  Although 
Na2lr03  is  not  a  quantum  spin  liquid,  the  fact  that  the 
ordering  temperature  =  15  K  is  considerably  smaller 
than  both  the  Curie-Weiss  temperature  r0  =  — 125  K  [7] 
and  the  spin  wave  energy  ~  5  meV  [10]  implies 
that  the  degree  of  frustration  is  still  quite  strong  and  the 
Kitaev  term  should  dominate  the  low-energy  physics. 
Nevertheless,  despite  intensive  research  performed  on 
Na2lr03  thus  far,  to  the  best  of  our  knowledge,  there 
has  been  no  evidence  of  the  spin-liquid-type  behavior  in 
this  material,  which  should  show  up  in  the  ^  T 
temperature  range  [14]. 

In  this  Eetter  we  report  our  results  on  ultrafast  studies  of 
photoexcitations  in  Na2lr03.  We  observe  a  change  in  their 
dynamics  across  Neel  temperature;  namely,  we  observe  a 
sharp  increase  in  the  binding  energy  of  the  excitons  that 
they  form  as  the  system  enters  the  ordered  phase.  We 
interpret  this  as  evidence  of  confinement-deconfinement 
transition  of  spin  and  charge  excitations  across  which  is 
a  hallmark  of  spin-liquid  physics  [15,16]. 

Time -resolved  experiments  were  performed  with  a  Ti: 
sapphire  oscillator  lasing  at  the  center  wavelength  of 
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795  nm  {hco  =  1.55  eV)  producing  pulses  of  60  fs  in 
duration.  The  repetition  rate  of  the  laser  was  reduced  to 
1 .6  MHz  with  an  external  pulse  picker  to  avoid  cumulative 
heating  effects  on  the  sample,  and  the  spot  size  of  a 
spatially  Gaussian  beam  was  set  to  60  ptm  EWHM. 
Single  crystals  of  Na2lr03  were  grown  using  a  self¬ 
flux  method  from  off-stoichiometric  quantities  of  Ir02 
and  Na2C03.  Similar  technical  details  were  described 
elsewhere  [17-19]. 

Data  were  obtained  with  a  standard  optical  pump-probe 
technique  [20,21]  where  a  single  “pump”  pulse  excites  the 
sample  and  the  resulting  dynamical  response  is  monitored 
by  the  normalized  change  in  the  reflectivity  AR{t)/R  of  a 
separate  “probe”  beam  as  a  function  of  time  delay  At 
between  the  pump  and  probe.  We  use  the  same  wavelength 
for  both  pump  and  probe  pulses.  Eor  phase-sensitive 
measurements  a  variation  of  the  pump-probe  technique 
called  the  “heterodyne  transient  grating”  (HTG)  method  is 
used  [22].  This  method,  unlike  the  pump-probe  technique, 
can  distinguish  between  components  of  AR  with  different 
physical  origins  [23].  Here  an  interference  of  two  pump 
beams  produces  a  spatially  modulated  excitation  pattern 
which  is  studied  by  a  probe  beam  diffracted  off  the  sample. 
The  diffracted  beam  is  then  heterodyned  with  an  additional 
beam  used  as  a  local  oscillator.  The  time  dependence  AR{t) 
of  a  multicomponent  system  response  obtained  in  this  way 
changes  as  a  function  of  the  phase  difference  (p  between  the 
probe  beam  and  the  local  oscillator,  whereas  that  of  a 
single-component  system  just  scales  proportionally  to 
cos(^)  (also  see  the  Supplemental  Material  [24]). 

Given  the  band  structure  of  Na2lr03  [26-28]  [Eig.  1(a)], 
the  absorption  of  a  pump  photon  with  energy  E  eV  causes 
electrons  to  transition  from  a  Teff  =  3/2  band  into  the  upper 
Hubbard  band,  which  is  the  only  accessible  level  for  the 
excited  electrons  for  this  photon  energy.  The  depleted 
valence  band  is  then  partially  refilled  through  relaxation 
of  the  photoexcited  electrons  and  partially  with  electrons 
from  the  lower  Hubbard  band  [Eig.  1(a)].  At  the  end  of  this 
relatively  fast  process  the  system  will  have  some  amount  of 
excitations  in  the  upper  Hubbard  band  (double  occupancies 
in  real  space  or  “doublons”)  and  an  equal  amount  of  holes 
in  the  lower  Hubbard  band.  The  excited  state  is  metastable 
as  the  optical  dipole  transitions  within  the  Hubbard  band 
are  prohibited  by  selection  rules  (AJ  =  0  transition). 
Moreover,  the  energy  of  the  magnons,  which  are  the 
relevant  excitations  c«5-10  meV  [10],  is  much 

less  than  the  Hubbard  gap  U  ~  350  meV,  which  is  the 
energy  which  needs  to  be  dissipated  during  the  doublon- 
hole  recombination  process,  making  the  lifetime  exponen¬ 
tially  large  in  U/e  [29].  This  observation  allows  us  to 
consider  holes  and  doublons  as  stable  quasiparticles  for  the 
time  scales  relevant  to  our  experiments  (~100  ps). 

Eigures  1(b)  and  1(c)  show  HTG  data  taken  at  a  pump 
fluence  of  ~9.5  ^i/crc?  for  various  values  of  cp  at  295  and 
25  K,  respectively.  At  each  temperature  the  shape  of  the 


Time(ps) 


Temperature  (K) 


FIG.  1  (color  online),  (a)  Sketch  of  the  band  stmcture  of 
Na2lr03  illustrating  the  relevant  processes  during  photoexcitation 
with  1.5  eV  light,  (b)  HTG  traces  for  selected  values  of  (p  (see 
text)  at  a  pump  fluence  of  9.5  jij/crc?  at  T  =  295  K  and  (c)  at 
T  =  25  K.  Note  the  multicomponent  behavior  of  the  HTG  traces 
as  a  function  of  (j)  and  the  relative  strengthening  of  the  initial 
spike  (labeled  with  an  arrow)  associated  with  single-particle 
excitations  (SE)  at  higher  temperature  (see  text),  (d)  Phase  {&)  of 
the  total  signal  at  t  k  50  ps  in  the  quasiequilibrium  state  as  a 
function  of  temperature.  Error  bars  represent  the  95%  confidence 
interval  (2  s.d.)  in  extracting  the  phase.  Solid  red  line  is  a  fit  to  the 
data  based  on  the  Boltzmann  distribution  of  SE  and  Hubbard 
exciton  (HE)  populations,  SE/HE  oc  exp(— A/kgT).  Inset: 
Phaser  diagram  representing  quasiequilibrium  AR  (black  phasor) 
due  to  SE  and  HE. 


differential  reflectivity  time  trace  AR{t)/R  with  time 
changes  as  cp  is  varied,  indicating  that  there  is  more  than 
one  component  in  the  system  response.  This  behavior  is  in 
agreement  with  earlier  resonant  inelastic  x-ray  scattering 
(RIXS)  studies  [26,30],  which  demonstrated  that  the  low- 
energy  excitations  of  Na2lr03  are  single-particle  (doublons 
and  holes)  excitations  (SE)  and  their  bound  state  is  known 
as  a  Hubbard  exciton  (HE)  [31].  Eigures  3(a)  and  3(b) 
demonstrate  that  the  component  featuring  the  fast  spike 
near  t  =  0  is  clearly  getting  stronger  with  increasing 
temperature,  allowing  us  to  identify  it  with  SE.  This  is 
similar  to  the  results  of  studies  of  photoexcited  Mott 
insulators  in  other  systems  [31]  where  the  AR/R  compo¬ 
nent  with  an  initial  fast  spike  was  shown  to  be  due  to  SE 
whereas  the  one  without  is  due  to  HE. 

It  should  be  noted  that  both  components  (SE  and  HE)  are 
long-lived  and  thus  they  both  contribute  to  the  composition 
of  the  total  signal  in  the  long-time  limit.  There,  a  quasi- 
thermal  equilibrium  is  estabhshed  between  SE  and  HE 
and  thus  their  population  ratio  should  be  proportional  to 
the  Boltzmann  factor  exp(— A/k^T),  where  A  is  the  HE 
binding  energy.  In  this  regime,  the  net  phase  6  of  the  signal 
response  reaches  a  constant  value  that  is  directly  related  to 
the  quasiequilibrium  population  ratio  of  SE  and  HE  (see 
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the  inset  of  Eig.  1(d)  and  the  Supplemental  Material  [24]). 
This  phase  d  is  plotted  in  Eig.  1(d)  as  a  function  of 
temperature  from  which  we  extract  A  «  4.6  ±  0.8  meV. 
This  value  is  within  the  hounds  set  hy  RIXS  measurements 
[26]  and  confirms  that  the  component  featuring  a  fast  spike 
at  t  =  0  is  indeed  due  to  SE.  Eigures  2(a)  and  2(h)  show  the 
HTG  data  taken  for  a  very  low  pump  fluence  of  value  of 
~30  nJ / cm^  at  temperatures  above  and  helow  T respec¬ 
tively.  Similar  to  the  higher  fluence  data,  the  low  fluence 
response  above  [Eig.  2(a)]  clearly  features  more  than 
one  component,  indicating  the  presence  of  both  SE  and  HE. 
On  the  other  hand,  the  low  fluence  response  below 
[Eig.  2(b)],  strikingly,  scales  proportionally  to  cos{c/)), 
implying  a  single-component  behavior,  which,  as  discussed 
above,  is  due  to  Hubbard  excitons. 

We  now  proceed  to  study  this  disappearance  of  SE  as  a 
function  of  temperature  by  performing  optical  pump-probe 
measurements  for  low  excitation  densities  to  minimize 
heating  effects.  Eigure  2(d)  shows  reflectivity  transients  for 
various  temperatures  below  for  two  different  pump 
fluences.  As  can  be  seen,  the  normalized  system  response 


FIG.  2  (color  online).  HTG  and  pump-probe  data  at  various 
temperatures  at  low  pump  fluence.  (a)  Representative  HTG  traces 
at  r  =  17  K  for  two  different  values  of  ^  at  a  pump  fluence  of 
30  nJ/cm^.  The  traces  exhibit  qualitatively  different  shapes 
indicating  the  presence  of  both  SE  and  HE.  Note  the  fast  spike 
around  ?  =  0  for  the  purple  curve.  This  is  signature  of  SE  (see 
Fig.  1  and  text),  (b)  HTG  traces  at  T  =  5  K  for  19  different  values 
of  (/i  at  a  pump  fluence  of  30  nj/cm^,  scaled  to  emphasize  the 
single  component  (HE)  nature  of  the  response.  Inset:  Unsealed 
HTG  traces,  (c)  Scaled  AR/R  traces  for  temperamres  above  Tn 
{T  =  17,  25,  and  35  K)  at  different  fluences:  100  nj/cm^  (filled 
markers)  and  50  nj/cm^  (open  markers).  Note  the  strong  temper¬ 
ature  dependence  and  the  lack  of  fluence  dependence  in  this  limit, 
(d)  AR/R  traces  for  temperatures  below  (T  =  5,  7,  10,  and 
12  K),  scaled  to  emphasize  the  universal  behavior  of  transient 
traces.  Upper  curve,  50  nj/cm^;  lower  curve,  30  nj/cm^.  Curves 
at  different  fluence  values  are  shifted  for  better  clarity. 


in  this  regime  is  independent  of  both  temperature  and  pump 
fluence,  demonstrating  that  the  single-component  behavior 
observed  at  T  =  5  K  persists  up  to  T^.  This  indicates  that 
SEs  are  suppressed  throughout  the  ordered  phase.  On 
the  contrary,  reflectivity  transients  above  [Eig.  2(c)] 
strongly  depend  on  temperature,  which  combined  with  the 
HTG  data  (Eig.  1)  indicate  the  formation  and  strengthening 
of  the  component  due  to  SE.  Moreover,  above  T^,  the 
normalized  transients  at  each  temperature  are  independent 
of  pump  fluence  [Eig.  2(c)],  demonstrating  that  the  relative 
composition  of  the  signal  (ratio  between  SE  and  HE 
populations)  is  constant  as  a  function  of  pump  fluence 
in  this  low  excitation  density  regime.  This  sudden  dis¬ 
appearance  of  SE  at  implies  a  sharp  increase  in  the  HE 
binding  energy. 

In  general,  an  increase  in  binding  energy  of  an  exciton 
can  be  either  due  to  an  enhancement  of  the  attracting 
potential  or  due  to  “slowing  down”  of  the  overall  dynamics, 
e.g.,  by  increasing  the  effective  mass.  The  effective  mass 
of  a  single  hole  in  a  Mott  insulator  is  indeed  enhanced 
due  to  emission  of  magnons,  but  this  happens  in  both 
the  antiferromagnetic  and  disordered  phases  [32].  But 
more  importantly,  it  is  the  sheer  disparity  between  the 
bandwidth  of  the  single-particle  excitations  (>100  meV, 
e.g..  Ref.  [26])  and  the  Heisenberg  couphng  responsible  for 
antiferromagnetism  k  I  meY)  that  makes  the 

slowing-down  scenario  unlikely.  Therefore,  we  conclude 
that  the  increase  in  binding  energy  of  HE  is  a  result  of 
enhancement  in  attraction  between  SEs. 

Here  we  present  a  simple  intuitive  picture  of  the  effective 
attraction  mediated  by  antiferromagnetic  ordering  respon¬ 
sible  for  additional  binding  energy,  based  on  the  Kitaev- 
Heisenberg  model.  This  Hamiltonian  naturally  gives  rise  to 
the  zigzag  order  [9]: 

Hkh  =  +  JhYs.  -Sj.  ( 1 ) 

(ij)  {ij) 

The  main  term  here  is  the  strongly  frustrated  Kitaev  term 
(r@  =  —125  K).  The  Heisenberg  term  lifts  the  degeneracy 
and  the  system  “freezes”  into  an  ordered  state  below  T In 
the  zigzag  phase  every  spin  finds  its  “Kitaev  partner”  and 
antialigns  itself  with  the  spin  of  its  partner  in  the  direction 
determined  by  the  orientation  of  the  connecting  bond.  The 
much  smaller  Heisenberg  term  (T^  =  15.3  K)  tries  to 
minimize  its  energy  under  the  condition  that  every  spin 
has  a  Kitaev  partner.  As  can  be  seen  in  Eig.  3(b),  the  zigzag 
order  satisfies  this  condition  for  all  bonds  except  those  that 
connect  Kitaev  partners. 

This  situation  changes  drastically  when  spinless  defects 
such  as  doublons  or  holes  are  introduced  into  the  system. 
These  particles  can  be  thought  of  as  topological  in  the  same 
sense  as  the  excitations  in  the  Rokhsar-Kivelson  dimer 
model  [33].  The  spins  that  have  lost  their  Kitaev  partners 
reorganize  the  surrounding  spin  order  at  the  expense  of 
Heisenberg  energy  [Eig.  3(c)].  Reoriented  spins  form  a 
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FIG.  3  (color  online),  (a)  Schematic  view  of  an  excited  sample. 
Green  dots  correspond  to  doublons  while  the  white  dots  represent 
holes.  The  shaded  red  area  marks  the  region  where  the  spins 
are  affected  by  the  reconstruction,  i.e.,  the  “string”  [see  (c)]. 
(b)  Sketch  representing  the  zigzag  ordered  low  temperature 
ground  state.  Red  bonds  connect  Kitaev  partners  while  gray 
bonds  have  a  spin  configuration  that  minimizes  the  Heisenberg 
energy,  (c)  Simplified  representation  of  the  restructured  state  with 
singlet  defects  (hole  or  doublon).  The  Heisenberg  energy  along 
the  bonds  highlighted  by  purple  ovals  is  no  longer  minimal.  The 
string  composed  of  such  bonds  must  begin  and  end  on  a  defect, 
(d)  Illustration  of  a  phase  diagram  of  a  modified  Kitaev  model 
with  a  first-order  quantum  critical  point  and  a  phase  boundary 
between  confined  (CF)  and  deconfined  (DC). 

string  terminating  on  the  other  defect.  The  energy  cost  of 
this  configuration  is  proportional  to  the  number  of  broken 
Heisenberg  links,  which  in  turn  is  proportional  to  the  length 
of  the  string  connecting  the  two  defects.  This  prohibits  the 
long-range  separation  of  the  defects  which,  in  the  case  of 
low  excitation  density,  will  predominantly  be  of  opposite 
charges  (doublon-hole),  leading  to  an  enhanced  binding 
between  them  in  addition  to  Coulomb  attraction.  This  is 
similar  to  the  picture  of  the  quark  confinement  in  high- 
energy  physics  [34]:  separation  of  defects  produces  a  string 
of  perturbed  vacuum  between  them  with  an  energy  propor¬ 
tional  to  its  length.  However,  since  breaking  of  this  string  in 
our  case  produces  a  pair  of  electrically  neutral  unpaired 
“dangling”  spins  (which  are  confined),  the  binding  energy 
between  doublons  and  holes  is  limited  by  the  cost  of 
breaking  a  Kitaev  pair,  which  is  of  the  order  of  Kitaev 
coupling,  J K,  10  meV  in  the  case  of  Na2lr03. 

This  is  consistent  with  previous  theoretical  works  on  the 
Kitaev-Heisenberg  model  where  it  was  observed  that  for 
sufficiently  weak  perturbations  the  Kitaev  spin-liquid  state 
persists,  but  as  the  extra  term  gets  stronger,  the  system 
enters  an  ordered  state  [4,5,35].  Formulated  in  terms  of 
spinons,  such  a  transtition  corresponds  to  the  transition 
from  a  deconfined  state  (spin-liquid)  to  a  confined  state 


(antiferromagnet)  [35,36].  This  is  a  first-order  phase 
transition  [35,36],  and  no  quantum  critical  region  is 
expected  above  the  quantum  critical  point  X  =  [37], 

where  X  is  the  relative  strength  of  the  Heisenberg  term  in  the 
Hamiltonian.  Therefore,  confined  and  deconfined  phases 
are  separated  by  a  simple  boundary  and  the  confinement- 
deconfinement  transition  can  be  observed  not  only  by 
tuning  the  strength  of  the  perturbing  term  2  at  T  =  0  as  in 
Refs.  [35,36],  but  also  by  going  across  the  ordering 
temperature  for  a  fixed  X>  X^  [see  Fig.  3(d)].  In  the 
confined  phase,  all  fractional  excitations  such  as  holons, 
doublons,  and  spinons  are  bound  to  each  other  and, 
conversely,  can  move  independently  in  the  deconfined 
phase  [16]. 

The  fact  that  the  zigzag  order  is  not  a  trivial  antiferro¬ 
magnet  was  evidenced  in  a  recent  work  of  Manni  et  al.  [38], 
where  partial  substitution  of  Ir  atoms  with  nonmagnetic  Ti 
atoms  resulted  in  the  formation  of  a  spin  glass  state  at  low 
temperatures.  This  has  a  natural  explanation  within  our 
framework,  as  Ti  sites  can  be  treated  as  static  spinless 
defects.  At  sufficient  concentration  they  will  be  connected 
by  spin  strings  as  in  Fig.  2(c).  Since  there  are  many 
different  ways  to  connect  different  Ti  sites,  it  is  natural 
to  expect  a  spin  glass  state  at  low  temperatures. 

In  conclusion,  we  performed  an  optical  pump-probe 
study  of  Na2lr03,  a  material  proposed  to  be  a  realization  of 
the  Kitaev  model.  We  observed  that  photoinduced  charged 
excitations  display  drastically  different  behavior  below  and 
above  the  Neel  ordering  temperature.  Namely,  the  binding 
energy  of  the  excitons  that  these  particles  form  undergoes  a 
sharp  increase  upon  entering  the  ordered  phase.  Based  on 
earlier  theoretical  studies  on  doped  Mott  insulators,  we 
conjecture  that  this  is  due  to  an  effective  attraction  brought 
about  by  the  antiferromagnetic  order  rather  than  because  of 
an  increase  in  effective  mass  of  quasiparticles.  We  argue 
that  this  attraction  is  a  manifestation  of  confinement  of 
spin  excitations  anticipated  in  the  ordered  phase  of  the 
Kitaev-Heisenberg  model.  Therefore,  the  change  of  behav¬ 
ior  that  we  observe  at  the  Neel  temperature  is  due  to  the 
confinement-deconfinement  transition,  providing  evidence 
of  spin-liquid-type  physics  in  Na2lr03. 
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In  our  previous  work  which  was  published  in  PRL  in  January  2015,  we  have  established 
that  the  interactions  within  Na-213  give  rise  to  strongly  frustrated  Hamiltonian  consistent 
with  Kitaev-Heisenberg  model,  which  among  other  things  gives  rise  to  spin-liquid  type 
behavior  in  the  disordered  phase  above  Neel  temperature.  One  of  the  open  questions  left 
in  that  work  was  the  nature  of  the  slow  dynamics  of  the  signal  associated  with  the  ordered 
phase.  A  detailed  investigation  of  this  question  should  have  important  implications  for  the 
very  little  understood  problem  of  the  non-equilibrium  Mott  insulating  phase. 

An  implication  of  our  work,  which  might  give  a  hint  on  what  is  going  on  during  the  slow 
(tens-of-picoseconds-scale)  evolution  of  the  system,  was  that  in  the  presence  of  defects 
(either  permanent  defects  or  transient  ones,  such  as  slow  photo-excitations)  the  low- 
temperature  phase  of  Na-213  should  exhibit  spin-glass  behavior.  The  slow  process  that 
we  observe  might  therefore  be  the  relaxation  of  a  spin-glass.  The  reason  why  the 
timescales  we  observe  in  our  pump-probe  experiments  are  short  for  a  spin-glass,  are  then 
due  to  the  fact  that  the  disorder  due  to  photo-induced  carriers  is  not  a  perfectly  quenched 
disorder. 

In  order  to  investigate  the  nature  of  the  low  temperature  phase  of  a  frustrated  Mott 
insulator  Na-213  we  have  built  a  special  three  pulse  pump-probe  setup,  with  two  beams 
for  pumping  and  one  for  probing.  According  to  our  design  one  of  the  pump  beams  has  a 
wavelength  of  800nm  to  make  a  connection  with  previous  experiments.  As  before  this 
pulse  couples  most  efficiently  to  electronic  degrees  of  freedom.  On  the  other  hand  the 
other  pump  beam  is  deep  inside  the  infrared  spectrum  (2100nm).  This  beam  energy 
corresponds  to  forbidden  frequency  range  for  electronic  transitions  in  Na-213,  and 
therefore  mostly  couples  to  non-electronic  degrees  of  freedom  such  as  phonons  and  from 
the  point  of  view  of  electrons,  only  produces  perturbations  in  temperature.  With  this  setup 
we  could  do  an  ultrafast  analog  of  classical  experiments  of  M.  Lederman  et  al.  and  M. 
Ocio  et  al,  on  spin  glasses. 
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On  the  figure  above  we  show  our  preliminary  result.  The  horizontal  axis  on  both  plots 
(main  panel  and  inset)  corresponds  to  the  time  delay  between  the  2100nm  pump  pulse 
and  the  900nm  probe  pulse  (we  chose  this  wavelength  over  usual  800nm  to  avoid 
excessive  pump  scattering).  The  inset  panel  shows  the  full  measured  time  range  and  the 
main  panel  shows  the  zoom  onto  the  region  around  the  time  of  arrival  of  the  2100nm 
pulse.  Different  curves  correspond  to  different  times  of  arrival  of  the  800nm  pump 
relative  to  the  2100nm  one.  The  first  thing  to  notice  on  the  main  panel  is  that  after  the 
second  pump  pulse  (2100nm)  the  system  “remembers”  its  previous  history  for  a  while. 
This  is  consistent  with  our  expectations  that  2100nm  pulse  only  induces  sudden 
temperature  increase  without  affecting  electronic  distribution.  But  the  most  important 
message  of  the  figure  is  that  although  the  delay  between  800nm  and  2100nm  pulse  can  be 
as  large  as  40ps  (see  inset),  the  history  does  get  forgotten  completely  within  less  than 
lOps. 

If  our  picture  correct,  then  what  we  observe  in  our  regular  pump-probe  experiments  is  the 
relaxation  of  a  spin  glass.  Moreover  in  our  three-pulse  measurements  we  observe  the 
“forgetting”  effect  of  a  spin  glass  in  real  time.  The  comparison  of  the  characteristic  time 
of  “forgetting”  with  the  magnitudes  of  available  physical  parameters  of  Na-213  further 
support  this  picture.  Indeed,  quasi-classically,  the  low  energy  dynamics  of  the  Kitaev- 
Heisenberg  model  can  be  thought  of  as  fluctuating  “hard”  dimers  formed  by  strong 
Kitaev  term  weakly  interacting  among  each  other  through  Heisenberg  energy.  Therefore 
it  is  the  Heisenberg  energy  (Jj^  ~  T^eei  ~  Ips)  which  is  responsible  for  the  dimer  dynamics, 
and  ultimately  the  relaxation  times. 


